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Abstract 

The eigenvalues of the Hessian associated with random manifolds are constructed for the 
general case of R steps of replica symmetry breaking. For the Parisi limit R — > oo (continuum 
replica symmetry breaking) which is relevant for the manifold dimension D < 2, they are 
shown to be non negative. 



Resume 

Les valeurs propres de la Hessienne, associee avec une variete aleatoire, sont construites 
dans le cas general de R etapes de brisure de la symmetrie des repliques. Dans la limite de 
Parisi, R — > oo (brisure continue de la symmetrie des repliques) qui est pertinente pour la 
dimension de la variete D < 2, on montre qu'elles sont non negative. 



PACS numbers: 05.20, 05.50, 02.50 



1 Introduction 

The behaviour of fluctuating manifolds pinned by quenched random impurities encompasses 
a rich variety of physical situations, from the directed polymers in a random potential (with 
a manifold dimension D — 1) to interfacial situations (D — d— 1, d being the total dimension 
of the space). The directed polymer problem, for example, is itself related to surface growth, 
turbulence of the Burgers equation and spin glass. The literature concerning various aspects is 
wide and florishing, a good deal of it being quoted in some recent reviews^], ^). The approaches 
taken have been quite diverse and it is not the purpose of this paper to describe them. This 
paper is a first step in an approach to random manifolds that uses field theoretic techniques 
in the replica formalism. Field theoretic techniques have been used with success (in their 
functional renormalization group variety) by Balents and Fisherj| to describe the large N 
(N — d — D is the transverse dimension) for the 4 — e dimensional manifold. In the context of 
directed polymers Hwa and Fisher^], using more conventional field theory, have emphasized 
the role of large but rare fluctuations. Alternately a unified description of the general random 
manifold has been developed by Mezard and Parisi[^](MP), where, via the use of the replica 
technique, they have laid down the first step (mean field and self-consistent Hartree-Fock level) 
of a systematic field theory 1 /N expansion^]. This is the approach we want to pursue here and 
this paper may be considered as a first sequel to MP. 

In this paper we study the eigenvalues of the Hessian associated with the 1/N expansion 
as laid down by MP. We do it in a general form by keeping discrete the R steps of replica 
symmetry breaking, and we show that the MP continuous solution (R = oo, as is appropriate 
for D < 2) does not generate negative eigenvalues. 

The general discrete expression derived here allows to control carefully the limiting process 
to the MP continuous solution. It will be also of use in further discussions of the solutions for 
the case D > 2 (and in particular D — > 4) an unclear case to date. 

In section 2 we summarize the MP approach down to the quadratic fluctuation terms 
defining the Hessian (see also j?]]). Section 3 is devoted to the introduction of a discrete 
Fourier-like transform on the replica overlaps. In section 4 we briefly discuss the equation of 
state at mean field level. In section 5 we discuss and exhibit the replicon eigenvalues. Section 
6-8 are devoted to the tougher case of the longitudinal-anomalous eigenvalues. 



2 Summarizing MP approach 

We thus consider the MP random manifold model defined by 



The intricacies of the first loop correction were briefly considered in further work| 
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-W(t3) 



dx 



E 

.M=l 



+ 2V(x,lu(x)) + fjw 2 (x) 



(1) 



where d = N + D is the dimension of the space, TV the dimension of £5, i.e. of the tranverse 
space, and D the dimension of the vector x, the longitudinal directions. The quenched random 
potential V has a gaussian distribution of zero mean and correlation 



V(x,aj)V(x',uj') = -5 {l) >(x - x')Nf 



N 



(2) 



with 



f(y) = 



(8 + y) 1 - 



2(1-7) 

where the function is regularized at small arguments by the cut-off 8. 

Averaging over the random quenched variable V via replicas replaces (mj by 



(3) 



dx 




N 



(4) 

In studying (Q) MP have used first a variational approach (identical to the Hartree-Fock 
approximation), which becomes exact in the large N limit. Then they introduced auxiliary 
fields with the advantage that it generates correction to Hartree-Fock, namely 

(i) The zero loop term embodying the effect of quadratic fluctuation, i.e. the logarithm of 
the determinant of the Hessian matrix 



(«) 



corrections in a systematic way. 



In this work we are interested in (i) studying the Hessian and its eigenvalues around the 
mean field (Hartree-Fock or variational) solutions. We leave the study of (u) for a separate 
publication^]. 

Following MP we introduce the n(n + l)/2 component auxiliary fields 



r a b(x) = —LJa(x)LU b (x) 



(•>) 



and its associated Lagrange multipliers s a b(x). We need to compute the generating function 



2 Note that in ref. [8] the case 1— RSB, relevant for short ranged correlations, is extensively studied both for its 
stability and its 1/jV corrections (given in closed form). 
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a=l ^ a<b 

xexp|-^^ J dxs ab (x) \Nr ab (x) - Lu a {x)Lu b {x)]^ 



> cxp ^ — 77 / dx 



,/-=! 



now quadratic in w's, hence after integrating over uj a 



P 2 



dxy]Nf (r aa (x) + r bb (x) - 2r ab (x)) 



Z" = I Y[d[r ab ]d[s ab ]exp[Ng{r;s}] 

a<b 



Q{r;s} 



dx [sa6(a;)r £l t,(x) + /3/(r aa (x) + r 00 (x) - 2r a6 (x))] + 5[s] (7) 



expS[s] = J ]~[ d[ij a ] exp i -~ ^ J dx [jj a (x) ((-V 2 + p)S ab - s ab (x)) uj b (x)] > 
Assuming an uniform saddle point (for large N) with 



r a b(x) = p a b + 8r ab (x) 



s ab (x) = a ab + Ss ab (x) (8) 

one can expand around it. MP thus obtain 
(i) The mean field contibution Q^°'{r; s} 

(u) The stationarity condition (vanishing of linear terms) that defines p and a (identical to 
the variational answer^]) 

• Stationarity with respect to s a f,(x): 

3 Only in the large TV limit. If N is finite, in the variational answer f(y) is replaced by /({/) 
TJ ^ m daa N / 2 - 1 e- a f(2ay/N) with / -> / as N -+ oo. 
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p 

where 

{G- 1 ) ab = (p 2 + fi)8 ab - a ab (10) 
• Stationarity with respect to r ab (x): 



(Tab = 2/3/' (p oa + p bb - 2p ab ) (11) 

C aa — — CT a b (12) 

(iii) The Hessian matrix of the quadratic form in 8s ab (after integration of the quadratic 
terms in 8r ab which in particular says that Ss aa = ~^2 a ^ b s ab, i-e. 



Z» ~ exp [NG°{p; a}] j JJ d[«h o6 ] exp -± ^ ^ <5 Sai) (p)M a6;cd (p)<5 Scd (-p) (13) 



Kb a<b p 

c<d 



?ao:ca 

M^ d ( P ) = - , -y2{G ac ( q )+G bd ( q )-G ad ( q )-G bc (q)}x 

<tj (Paa + Pbb — tPab) ' 

x {G ac (p -q) + G bd (p -q)- G a d( P - q) - G bc (p - q)} (14) 



3 Discrete Fourier Transform 

We shall work in a formalism where the level R of replica symmetry breaking is kept discrete 
(R = replica symmetric, R = 1 one step of breaking, R = oo Parisi breaking, i.e. continuous 
overlaps). In the end, if necessary, we shall let R — > oo. 

For an observable o a/ 3, with a a p = a r for an overlap a n /3 = r, define the discrete Fourier 
transform as 



ij+i 

Ofc = 

r — k 

and its inverse 



^p r (a r -a r _i) (15) 



a r = ^ — (afc — afc+i) (16) 



fc=o 



G 



Here the p r 's are the size of the successive Parisi boxes, po = n, pn+i = 1. The above transform 
has the properties 



a r — a r -i = — (d r — a r +i) (17) 

Pr 

^^a a7 fo 7/ 3 = C a fi (18) 



a k b k = c fe (19) 



In particular 



(a- 1 )* = Y (20) 
All this is but a discrete version of the operations introduced in MP for the continuum (R — * 
oo). Note that in the above all quantities with indices outside the range (0, 1, . . . , R + 1) arc 
to be taken identically null. If the limit R — * oo is understood as 



a r — > a(u) 

a r — > a(u) (21) 
then the connexion with MP is, in their notations, 

a(u) — a — (a) — [a](u) 
[a](u) = - dta(t)+ua(u) (22) 



4 The equation of state 

With the above definitions we have 



V 



Gk( P ) = -5 , I r (24) 



P 2 + [I - CTfc 
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and the equation of state 



or = 2/3/' [ G R+^l) - G -(l)]j 



(25) 



= -^p t (ff t -(7 t _i) (26) 
The last equation, identical to (|l2]), also writes 

ffo = 53pt(ff t -ff t _x)=0 (27) 

r=0 

To solve for the equation of state one uses the transform that allows to write 

[G r (q) - G r -i(q)] = - [G r (q) - G r +l(q)] (28) 



Pt 



and 



r = \ [ G «+i(?) - G -(9)] = | \ £ ^-(» k -a k+1 )G k (q)G k+1 (q) + Gn+i(Q) 

q \k=r+l 



(29) 



The equation of state then writes, with 



a r = 2/3/ V) 



(30) 



i.e. 



1 r 1 

— 5 / , — [*fc - ^fc+i] 

k—u 



Here we have 



£ ^(ok-Ok+i)G k (q)Gk+i{q)+G R+1 {q) 



(31) 



G fc («) 



(32) 



3 2 + (i - a k 

where fj, is an infrared cut-off, which can be set to zero in appropriate situations. In the 
continuum limit 
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-oh -> \<y\(x) 



Pk — > a; 

(in the Parisi gauge) easily restituting the MP solution as given for the case of noise with long 
range correlations. 

For later use, let us take a one step difference in the equation of state 

Or - ffr-1 = — (ov — <5>+i) = 2/3 \f{y r ) - f'iUr-l)] (33) 

Pr L J 

with 

Vr-l ~Vr=lj2 _ G r _i(«)] = | £ ~ IA(«) " &+l(«)] (34) 



a 9 



i.e. 



2 [ Or — r+1 



J2Gr(q)Gr + l(<l) (35) 
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For small y, — i — y r (infinitesimal in the continuum) one gets 

l = -4j2Gr(q)Gr+i(q)f"M (36) 

9 

a formula used below. 

5 Hessian Eigenvalues: the replicon sector 

The Hessian ([jj) written in terms of overlaps becomes 

M ab ^ d (p) = M:: s v (p) (37) 

with 

r — af~]b , s — end 
u = max(a n c ; a D d) (38) 
v = max(b n c ; bC\d) 
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In ultrametric space, only three of the four overlaps are distinct. In the replicon subspace 
r — s, and we have two independent overlaps u, v > r + 1. In the longitudinal anomalous 
sector, if r 7^ s, one can keep max(u, v) to parametrize M. 

For the replicon sector, and on general ground [^, |l^, pj| , the eigenvalues are given by the 
double discrete Fourier transform 

R+l R+l 

u — k v= k 

for k,l > r + 1. 

Using the Hessian expression 

MTZ = J^+} 5 ^ R + 1 _ [G u (q) + G v (q) - 2G r (q)] [G„(p -q) + G v (p - q) - 2G r ( P - q)] 

(40) 



' ' 2/;' 



Hence, for the double transform 

K(r; k -. l ) = -JJ77-J2 [ G ^)G*(P - q) + G^Gkip - q)] (41) 



1 



and using (|3f]), valid in particular in the continuum, 

A p (r; k,l) = {2Gr(l)Gr+i(q) - [G k (q)Gi(p - q) + 6i(q)6 h (p - q)] } 

l,k > r + l (42) 

The propagators G r (q) are decreasing functions of |g| and of r (since — a r is an increasing 
function). Hence the most dangerous eigenvalue is 

A p (r; r + 1, r + 1) = 2 ^ {G r (g)G r+1 (g) - G r+1 (g)G r+1 (p - g)} (43) 
or in the continuum limit 



X p (x; x,x)-2 (2n)D Jdq D ^ 2+ ^_ ^ {p _ g)2 + fi _ &{x) J q2 + ^ _ a{x) (44) 

For D < 4, one can set the ultraviolet cut-off to infinity and rewrite the dangerous replicon 
eigenvalues, in the continuum limit, as 
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X p (x;x,x) = J2[G4l)-G x (p-q)] 2 (45) 
where the marginal stability is obviously displayed^]. 

6 Hessian eigenvalues: the longitudinal anomalous 
sector 

The analysis of longitudinal anomalous sector is notoriously more difficult. It has been con- 
siderably simplified by the introduction [12] hd, hll] of kernels, i.e. by going to a block-diagonal 



representation of the ultrametric matrices. Take in ( |37| ) a component of the longitudinal 
anomalous sector M£i£, r / s, which, given that only three overlaps can be distinct, can be 
written^] 

M t r ' s t = max(w; v) 

To M t r,s is associated a kernel (in fact a Fourier transform) K T k ' a , where k = 0, 1, 2, . . . , R + 
1, k = corresponding to the longitudinal sector. Likewise for the inverse Af _1 of the 
mass operator matrix we have a kernel F^' a . The relationship between them is the Dyson's 
equation Jlc|, pd[ 



R+l 

K" = - \ K^ffrK' 3 (46) 



£ 

t=0 

F[ is = -A fc (r)F- s A fc ( S ) (47) 



Here the Afc(r) are special replicon eigenvalues defined by 



X(r;k,r+1) k>r + l 
A,Ar)={ (48) 
X(r;r + l,r + 1) k < r + 1 



with A as in (B9|) and 



4 It is worth mentioning at this point that the limit D — > 4 is delicate to take. Indeed solving for the equation of 
state (30-32), one finds that the breakpoint x\ jjcx](a;) = for x > xij crosses its boundary value x\ = 1, at some 
intermediate dimension Do (7), 2 < Do < 4. This delicate point is left out for further consideration. See also Jlil 

5 Except for the components M^'v that also incorporate a replicon contribution. 
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A k (t) 



|(Pt — Pt+i) t < k — 1 

i(pt-2p t+ i) t = k-l (49) 



(pt-Pt+i) t>k-l 
The kernel K k ' 3 itself is a transformed of the mass operator matrix M t r,s . It writes as a 
generalized discrete Fourier transform 

(r s R+l \ 

E + 2 E + 4 E p ( M * r;s - M *- S i) ( 5 °) 
t=fc t=r+l i=s+l/ 

where we have displayed the case k < r < s (for different orderings the summation , 
is only kept in the allowed summands). Knowing K r h ,a , one then computes F£' s via (^, 
and (M _1 )(' s via an inverse transform. That is if one is able to solve the integral equation 
(|4(j) which is only possible with especially simple kernels K k ' 3 . This turned out to be the case 
when studying the bare propagator of the standard spin glass jlo| where 

\K^ = A k {min{r ]S )) (51) 

The same functional form also appears for the spin glass model introduced by Th.Nieuwenhuizen 
For the system studied here, let us construct the kernel K r k ' 3 with the matrix M^ 3 as given 
by the Hessian (jl4|) manifestly symmetric in r, s. Let us choose r < s and we then find 

R+l 

-K r k ' 3 (p) = -J2 E Pt{[G t (q)G t (p-q)-G t - 1 (q)G t -i(p-q)]- 

q t— max(s+l,fc) 



-G.{q)[G t (p-q)-Gt-i(p-q)]- [G t (q) - G t -i(q)]G 3 (p - q)} (52) 

which also rewrite as 

B k (s) r<s<k-l 
= <! (53) 



Bs+i(s) s>k-l 



where 



Bk(s) = - { ( G (9)5^ - «)) k - Gs(q)G k (p -q)- G k {q)G s {p - ?)} (54) 

9 

The important point is that the kernel only depends on max(r,s) 
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i^ 3 = B fc (max(r; S )) (55) 



thus rendering soluble the integral equation (ji^) . 



7 Solution of the Dyson equation: the propagator 
kernel 

In the case of \K r k ' a — y4ft(min(r; s)) the solution of the Dyson equation (j^j has been given 
(in the continuum limit) in We wish to extend it here to the case \K^' a = _Bfe(max(r; s)). 
In this section we work for a given value of k, i.e. we remain in a sub-block of the longitudinal 
anomalous sector (k = for the longitudinal). 

The propagator kernel of eqs ([l6| ^) is obtained as 

1 F r-.s = <t>j (min(r;s))^(max(r;s)) 

4 k A k (r)A k A k (s) { ' 
where is a Wronskian and 0^ the regular and irregular solution of a discrete Sturm- 

Liouville equation and the special replicon eigenvalues of (|48|). Namely we have for the 
Wronskian 

Afc = v^d) (57) 

with 

T> s f(s) = fs+i - f s (58) 

The functions <f> are given by 



with -Bfc(s) as in (^). Note that, contrary to the solution obtained for (|5l|), where Cfe = 1, 
here, </>^ and </>^" remain coupled via Ck 

i=0 

In terms of cpt the Wronskian can also be written as 
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1 + 



2>o 



(61) 



8 Stability discussion in the LA sector 

Given the explicit expressions (^) of Bk(s) in terms of G t (or Gi) itself a positive definite 
function, decreasing with increasing t 



one can easily show that 



V t Gt 



1 o"t+i — <5"t+2 



(62) 



(i) B(r) < 



(63) 



(m) © r B(r) = B(r + 1) - B(r) > 



(64) 



Under these conditions we can now give a lower bound for the longitudinal anomalous eigen- 
values. We can write out the determinant of the eigenvalues restricted to the sub-block k 



det (Af£> - Xl) = Y[ [Afe(r) - A] 



A fc ft) 
Afc(t) — A 



(65) 



which is directly related to the Wronskian (apart from the Ck factor). 
Note that, from equation (p9), one also obtains 



V s <f>+(s) = -~2V a B k (s)Y^ 



(66) 



With the boundary value c/>t(0) = 2 and T> s B(s) > 0, one infers that, in case of A < Ak(t) for 
all t, 



<t>t(s) > 
Vs<t>t{s) > 



(67) 
(68) 



Consider now in equation (|6^) a small value of A such that A < A*, (t) for all t . For that value 
A the determinant of ( p5| ) cannot vanish, it is always positive on account of ( |63| , |67j). Hence 
the spectrum of the longitudinal anomalous sector is certainly bounded by the bottom values 
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of the replicon spectrum A p= o(r;r + 1, r + 1)- This establishes stability, a marginal stability 
as in the continuum limit where A p= o(a;; x,x) = 0. 



9 Conclusions 

We thus have set up for random manifolds the formalism to discuss the stability of mean field 
solutions. For the well established continuous (R = oo) solution, adequate for D < 2, we are 
able to prove its marginal stability. The formalism used here allows a discussion of any type 
(discrete, continuous or mixed) of mean field solutions. 

It is worth noticing that, when discussing the stability, the only place where properties of 
the noise correlation function may play a direct role (outside of the fully continuous R = oo 
solutions) is in ( [33| - |36| ). There the equation of state is used to find the alternative analytical 
form for /", that cast in evidence the non-negative character of the replicon eigenvalues. 
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